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Three tasks were defined in the SOW; (1) Spectral Fitting of Current Data, (2) Shock Layer 
Modeling Using DSMC Code, and (3) Production of Nitrogen Oxides and Excited Nitrogen in Flow 
Fields. All three tasks were carried forward satisfactorily over the given time period with 
some alteration from the initial set of subtasks. Some changes due to the nature of the research, 
approved by the technical monitor, provided results beyond the scope of the original proposal. 
An abstract was completed, resulting from efforts on task 3, and the full paper, scheduled for 
presentation at an AIAA meeting in January, is near its final form. 


Task 1 - Spectral Fitting of Current Data (320/110 hrs) 

Proposed . It was proposed that analysis of N 2 emission from 8.72 MJ/kg shock layer be made 
at 2.54, 1 .91 , and 1 .27 cm positions and that vibrational state distributions, temperatures, 
and relative electronic state populations be obtained from data sets. It was further proposed to 
review other recorded arc-jet N 2 and air spectral data and to study NO emission characteristics 
characteristics in data. Spectra of selected data sets was to be analyzed for an abstract and 
technical report or journal paper. 

Accomplished . It was possible to adjust analysis techniques such that relative populations of 
two electronic states, the state of the neutral molecule and the A 2 II state of the molecular 
nitrogen ion could be determined for four positions within the 8.72 MJ/kg shock layer. Results 
show a change in population distributions and overall relative population between the two 
electronic states as one moves from the front of the shock to the boundary layer. Also, 
preliminary data has been obtained on relative populations between the first two excited 
electronic states of the ion , A 2 II and B 2 L. Such population ratios, an example of which is 
given in Appendix A, will provide information on electronic state distributions or electronic 
temperatures. Use was made of the PSI code over a range from 350 nm to 800 nm was made for 
these results. 

Evaluation of accuracy in the overall analysis was thought necessary and additional time 
has been used in evaluating errors. Since vibrational temperature analysis techniques require 
rotational temperature input, rotational temperature was again investigated using a different 
method of evaluation. The method was based on that used by Leger, et al. in a nitrogen plasma jet 
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study reported from France in 1990. Rotational temperature analysis is made using the NEQAIR 
code and new results show differences on the order of 300 K from previous measurements. This 
represents differences of only 4-8%, about equivalent to the accuracy of analysis. Additional 
information may be found in Appendix B. 


Other recorded arc-jet N 2 and air spectral data were reviewed and a technical report covering 
this data initiated. Part of that information is displayed in Appendix C. There are six air data 
sets, all taken near the front of the shock. Of the three N 2 data sets, two sets were taken at the 
improved spatial resolution and only only one of those, that from which much of the previous 
analysis has come, covers four points within the shock layer. 

Recommended . For improvement in the accuracy of analysis where long wavelength spans are 
necessary a change in the NEQAIR code from 20,000 to 40,000 points is recommended. But 
Park's updated code is likely set at 40,000 and likely will run on the cray without the 
modification of our current code. The A 2 II state it should be added to that code. Excitation as 
well as radiation data will be required. For the type of measurements desired the PSI code 
requires improved accuracy in radiation data for states above V = 12. 

Task 2 - Shock Layer Modeling Using DSMC Code (173/60 hrs) 

Proposed . A review of operational procedures of the then current JSC DSMC code was proposed. 
Information on other appropriate codes and modifications, including ionization, were to be made, 
and a determination of the applicability of codes reviewed to task requirement and comparison of 
computational capabilities and input variables was proposed. 

Accomplished . Two versions of the JSC DSMC code were found and their computational 
procedures were reviewed. An additional variation of the Bird DSMC code format was 
subsequently obtained from Sandia. The JSC code operated on the VAX had been modified for time 
sequence file output for particular plot considerations. Concentration was eventually made on 
the other JSC code operation from the cray. Studies were first made using simple boundary 
parameters for a spherical body as test article. 

New region and cell boundary parameters were identified by Scott for a flow investigation of a 
mass spectrometer probe system. The configuration is given in Appendix D. It was basically 
shown that the computational procedures could be used for this application, a relatively complex 
configuration, but some detail changes and additions to the code, such as wall catalysis, are being 
reviewed in order to make the results more useful. A plot routine placed on the Iris workstation 
was determined most ideal for the type of data reduction required, but no significant plots were 
obtained for inclusion in this report. 


n a commended . Acquisition and modification of a code which includes ionization and radiation 
from Langley or a developmental effort is recommended. 


Task 3 - Production of Nitrogen Oxides and Excited Nitrogen in Flow Fields 

(976 hrs) 

Proposed . It was proposed that a review be made of computational procedures used in CFD codes 
of Li and other codes on JSC computers and that the literature related to task objectives be 
reviewed. Analysis of problems associated with integration of specific chemical kinetics 
applicable to task into CFD codes was to be made. 

Accomplished . Computational procedures and literature related to the task objectives were 
reviewed. Results of this initial task led to a greater emphasis being placed in the analysis upon 
computation of basic model parameters required for CFD input, the excitation and relaxation of 
the nitrogen molecule. 


The diffusion equation, derived from the master equation for the electron-impact vibrational 
relaxation based on the diffusion approximation, was numerically solved for evaluation of time 
variations of vibrational number densities of each energy level and of the total vibrational 
energy. Subsequently, an approximate rate equation and a corresponding electron-impact 
vibrational relaxation time from the excited states was derived, compatible with the system of 
flow conservation equations. This represented the first such derivation considering multiple- 
quantum transitions from excited states to higher/lower states for the electronic ground state of 
the nitrogen molecule. 

The empirical curve-fit formulae for the improved relaxation time, which will be convenient 
for actual CFD flow calculations, were obtained using two quadratics in the logarithm of the 
electron temperature. The solution of the master equation for the e-V process shows a non- 
Boltzmann population distribution at the earlier stage of relaxation, which in turn suppresses 
the equilibration process. These results will be presented at the AIAA 30th Aerospace Sciences 
Meeting at Reno in January, 1992, as a paper titled Electron-Impact Vibrational Relaxation IP 
Hinh-Temperature Nitrogen . The abstract is included as Appendix E. 

Rpr.nmmanded . Further development of the approximate rate equation, which accounts for the 
non-Boltzmann nature of the population distribution is recommended. Analysis should be made 
for the rotational excitation process and the electronic excitation distribution of N 2 molecules at 
very high temperatures. 


APPENDIX A 


Preliminary Population Distributions at 3.18 cm 
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Preliminary Population Distributions at 2.54 cm 
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The population data shows first zeros populations for some vibrational levels indicating 
problems and unreliability of the numbers, those from the N 2 and N 2 + (A)-states 
particularly. The A-states of neutral molecule and ion produce radiation within the 
same wavelength range, so the problems with the fit are understood, requiring 
refinement in the process and additional vibrational data. The nonzero numbers, 
however, although inaccurate, represent reasonable relative populations for the 
vibrational states. 

A rough estimation of relative populations of the electronic states may be made, for 
example, by looking at the vibrational state populations for v' = 2. That would yield 

N 2 (A)/N 2 + (A) = 0.7 at 3.18 cm 

N 2 (A)/N 2 + (A) = 1.1 at 2.54 cm 


N 2 + (A)/N 2 + (B) = 807 at 3.18 cm 
N 2 + (A)/N 2 + (B) = 803 at 2.54 cm 

The latter two ratios would indicate little change in the electronic temperature of N 2 + 
between the two positions, but the first two indicate an increase in the relative 
population of the neutral molecule. 

Additional analysis as well as greater accuracy is necessary before definitive arguments 
will be made. The additional data and analysis will be presented in the completed paper. 


APPENDIX B 


One problem in determining rotational temperature is the overlaping of the P and R 
branches in the first negative system of the ion where the measurement is usually made. 
In such cases the overlapped lines define global lines from which it is possible to define a 
rotational temperature. Following the method used by Leger, the line intensity 
distribution as a function of temperature, t, would be 

8(X,t) = £p(X) + 6 r (X) 

a sum of integrated intensities of rotational lines of the P and R branches falling within 
an unresolved wavelength range with a peak at wavelength x. With much of the bandhead 
unresolved and consisting primarily of lines from the P branch, a bandhead intensity is 
defined by 


^head ^ 

where K is the total angular momentum quantum number apart from spin. Using the 
ratio, 



I(.) = 

^head 

for both computed, in our case by NEQAIR, and experimental spectra, a function 

E(t) = tL [ ( 1 < ex p) • i ( cal) ) 2 p 

i = 1 

using a square of differences is defined. Here N is the number of global lines, six in our 
case. Executing the least squares technique, with NEQAIR providing the calcuated 
intensities as a function of equilibrium temperature, the measured rotational 
temperature is that which minimizes E(t). 



Previous measurements utilized a different spectral matching and ratio technique. Those 
measurements have been repeated using the above technique and a comparison between 
results is given below. 


Rotational Temperature - two methods 
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The lower line in the graph represents the older measurements. Subsequent analysis of 
the new measurements show a sensitivity to careful matching of computed spectral 
characteristics to experimental spectra. 

A view of the spectral lines, as given in the series of data which follow, shows that the 
peak intensities of the rotational lines do not follow a normat temperature distribution. 
Consequently one cannot guess at a temperature by looking at the distribution. Although 
it appears to be a fluctuation phenomena which causes some error, temperatures can be 
inferred. 
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APPENDIX C 


ARC JET DATA SETS 


LABEL 
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GRATING 
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cm 
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7412 

air 

1100 

0.08 



600 

D 

7408 

air 

600 

0.04 


3.5 

600 
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800 

0.08 


3.5 

600 
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8333 
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800 

0.08 


3.5 
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8333 
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0.08 
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0.12 
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0.12 
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The most significant analysis has been made using the N- senes, where data was collected 
from four points within an N2 shock layer. 


The following spectra were taken at a low flow rate, 0.04 Ibs/s, with an arc current of 
600 A. These spectra correspond to that labeled D above. 
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The spectrum between 300 and 800 nm shows a single strong atomic line at 760 nm. A 
closer look at the spectrum below 700 nm is shown below. 
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While the first positive N2 bands appear clear for analysis, the first negative N2 + bands 
are obscured by NO radiation. To analize the N2 rad'ation bands for an accurate 
vibrational temperature or population distribution, a rotational temperature is needed 
(estimates can be made, however, through spectral matching in the analysis of 
vibration). Although analysis of the shortwavelength radiation, shown in greater detail 
in the following three spectral plots, appears exceptionally difficult, improved tools 
may allow significant information to be obtained. 
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Radiation overlaping the N2+ bandheads in the spectra above can be seen. 



The bandheads at 800 A, 0.08 Ib/s. labeled E, appears clearer as the 2nd positive 
system of N2 replaces the NO radiation. 
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Differences can be seen in the spectra labeled C with current increased to 1100 A at 
0.08 Ib/s. 


TITLE 1 
TITLE! 
TITLES 



c 


TITLE! 

TITLE! 

TITLE3 



c 



RELATIVE NTEN9TY *10® RELATIVE NTENSTY 


TITLE 1 
TITLE! 
T1TLE3 



THLEl 

TfTLE2 

TnLE3 







TITLE 1 
TITLE 2 
TrTLEJ 



WAVELENGTH, (nm) 

c 


TrTLEl 

TFTLE2 

TTTLE3 



c 


Information may be obtained also from the atomic lines. 




The air spectra labeled K taken with the flow rate increased to 0.12 Ibs/s are shown 
below. 
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APPENDIX E 


On the following pages is a draft of the paper on vibrational relaxation to be presented at 
the AIAA Aerospace Sciences meeting in Reno in January, 1992. 


Draft Prepared for AIAA 30th Aerospace Sciences Meeting 
January 6-9, 1992, Reno, Nevada 


ELECTRON-IMPACT VIBRATIONAL RELAXATION in 
HIGH-TEMPERATURE NITROGEN 

Jong-Hun Lee 1 

BSA Services, Houston, Texas 


Abstract 

Vibrational relaxation process of Nj molecules by elecrton-impact is exam- 
ined for the future planetary entry environments. Multiple-quantum transi- 
tions from excited states to higher/lower states are considered for the electronic 
ground state of the nitrogen molecule N:(X l E+). Vibrational excitation and 
de-excitation rate coefficients obtained by computational quantum chemistry 
are incorporated into the "diffusion model” to evaluate the time variations of 
vibrational number densities of each energy state and total vibrational energy. 
Results show a non-Boltzmann distribution of number densities at the earlier 
stage of relaxation, which in turn suppresses the equilibration process. An ap- 
proximate rate equation and a corresponding relaxation time from the excited 
states, compatible with the system of flow conservation equations, are derived 
for the first time. The relaxtion time from the excited states indicates the weak 
dependency of the vibrational temperature, and is shorter than the previously 
obtained relaxation time in which only excitation from the ground state was 
considered. The improved rate equation and the relaxation time, suited for 
the numerical simulation of the highly ionized planetary entry flowfields, are 
suggested. 


Nomenclature 


E t = total energy 

E v = vibrational energy 

E* = equilibrium vibratinal energy at electron temperature, T e 
f = correction factor 

h = Planck’s constant 

i Research Scientist. Mailing Address: NASA Johnson Space Center, Houston, TX 77058 
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J = rotational quantum number 

k VfV t = rate coefficients from vibrational state v to v' 

Ko = equivalent diffusivity 

M 2 = second moment of vibrational transitions 

m = highest vibrational state 

in e = electron mass 

Ti e = number density of electron 

TiM = number density of molecule 

rit, = number density of vth excited state 

Ti v » = number density of v'th excited state 

Pg — pressure of electron 

3 = arbitary parameter in rate equation 

t = time 

T = translational temperature of heavy particle 

Tg = electron temperature 

T v = vibrational temperature 

v } v 1 = vibrational quantum number 

6 V = characteristic temperature for vibration 

k = Boltzmann constant 

if = frequency of oscillator 

p v = normalized number density of vth state 

r e = vibrational relaxation time for e-V process 


Introduction 

Recently, there is an intense worldwide interest on the prosepct of both un- 
manned and manned exploration of the Moon and the Planet Mars. In order to 
reduce the weights of vehicles for such future missions, an aerobraking maneuver 
in'the atmosphere is highly desirable. The aerobraking maneuver takes place 
at the upper atmosphere at very high velocity over extended periods of time. 
As a result, the flow around such vehicles is expected to be both thermally and 
chemically in a nonequilibrium state. 

While the estimated entry velocity of the aeroassisted orbital transfer vehicle 
(AOTV) is about or less than 10 km/s, the Lunar-return entry velocity into the 
Earth’s atmosphere is up to about 11 km/s. The proposed Mars-return entry 
velocity in manned mission is between 11.3 and 14 km/s 1 . Unlike the flow 
around the AOTV, the flow surrounding the future aeroassited planetary entry 
vehicle (APEV) will be highly ionized because of higher entry velocities 3 . 

Due to its light mass, electron-molecule collisions are much more frequent 
than heavy particle-molecule collisions. Thus electrons are very efficient colli- 
sions partners in causing vibrational energy transitions in molecules. In electron- 
vibration (e-V) energy exchange processes, the Landau-Teller conditions, in 
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which only the transitions between adjacent levels are allowed t are violated 3 . 
Therefore, the multiple-level vibrational energy transitions must be considered 
in the e-V process. 

For atmospheric molecules, the low-energy electron-impact vibrational exci- 
tation process is usually dominated by resonant mechanism 4 . That is, an inci- 
dent electron is temporary captured in the neighborhood of a target molecule, 
thus forming a compound state. This resonance enhances the vibrational exci- 
tation cross sections by orders of magnitude. Lee 5 theoretically obtained the vi- 
brational excitation cross sections and rate coefficients for the electronic ground 
state of the nitrogen molecule with the boomerang model of Dube 

and Herzenberg 6 . Although this method is in principle a semi-empirical one, 
the obtained cross sections agree reasonably good with the experimental data 
of Schulz 7 . 

Lee 5 obtained the analytical solution of the time variation of vibrational 
energy based on the ’’diffusion model 8 ” at high-temperature regime. Subse- 
quently, an approximate rate equation and a corresponding relaxation time, 
compatible with the system of flow conservation equations, were developed. 
Although the obtained relaxation time agrees substantially with experimental 
result as pointed out by Hansen 9 , only the multiple-level excitation from the vi- 
brational ground state is considered in the evaluation. This assumption is valid 
only for the flow in the vicinity of the shock wave around an vehicle, where 
most molecules are populated in the vibrational ground state (v = 0) at t = 0. 
However, if the flow leaves the’shock wave, the contributions of the vibrationally 
excited states (t> > 0) become important and the effects of excitation and de- 
excitation from these v > 0 states must be considered. This is especially true 
for the flow around APEV in which ionization level is much higher than that of 
AOTV. 

The purpose of present paper is to solve numerically the master equation for 
vibrational transitions from the excited states, and to develop an appropriate 
rate equation and a relaxation time which approximately represent the rigorous 
solution of the master equation. The obtained rate equation and the vibrational 
time for e-V process will be suitable for the numerical simulation of APEV 
flowfield, since they include the effects of multiple excitation/de-excitation from 
and to non-ground states. 


Vibrational Rate Process 
Master Equation for Vibrational Transitions 

Assuming that the vibrational oscillators interchange their energies with an 
electron heat bath having a constant temperature T ei the rate of change of the 
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number density n v without dissociation is given by the master equation 




dt 


v'=0 


( 1 ) 


where n v is the number density of vth vibrationaily excited state, k V)V > the rate 
coefficient for the transition from state t; to ti', m the highest vibrational state 
considered, and n e the number density of electron. The principle of detailed 
balancing provides the condition 

k v f tV 7l v *E == ( 2 ) 


where n v s is the equilibrium population in state r. Using the relation (2) 
and the normalized number density p v = n v /n v Ei one obtains the following 
normalized master equation 


(3) 

dt 

At the high tempeature regime where the kinetic energy kT c is larger than the 
vibrational energy gap, collisional excitation and de-excitation of the vibrational 
states occur almost in accordance with the classical mechanics 8 . According to 
the classical analogy, the vibrational energy states are continuously distributed 
over energy state v } and the summations in Eq. (3) can be replaced by integra- 
tions 

-%= f K,,.{P*-P*)dv' (4) 

Tl e OX Jv '— 0 

If Jk V)1f ' in Eq. (4) is assumed to be large only for |t/ — u| < 1, i.e., the transitions 
occur mostly between neighboring states, then the the integration limits can be 
changed from -oo to oo 

-i ?£. = n k v ,..(p v . ~ Pv )dv' (5) 

7l< dt oo 


The right-hand side of Eq. (5) is expanded in powers of d n p v /dv n , and the terms 
only to second order are retained , then, one obtains the following diffusion type 
equation 8,10 : 

= ( 6 ) 

n t dt dv ^ dv 


where 



(7) 


is the second moment of vibrational energy transfer between rth state and i/th 
state. When k v>v > is given for discrete vibrational states v and t/, the moment 
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M can be determined by summing over the the discrete states 

M(v) = -v) 2 k v y (8) 

The resulting equation (6) is a diffusion equation in the vibrational space 
v. This equation, together with the requirement that the sum of all n^’s equals 
the given number density of molecules 71 ^ 

TO 

( 9 ) 

»o Jo 

gives the number density n vy which is a function of time and energy state v. 


Solution of Master Equation 

In the following analysis, several assumptions are made. First, the harmonic 
oscillator model is applied to describe the molecular vibration. Second, the 
effects of dissociation are neglected. These assumptions are valid rigorously 
only for the lower vibrational energy states. High vibrational states deviate 
significantly from the harmonic oscillator assumption due to the anharmonic 
vibrational potential and the very small unequal energy gaps between the vibra- 
tional states 10 . Furthermore, dissociation of molecules depletes the population 
of upper vibrational states. In gas mixtures with electrons such as the case of 
planetary entry environments, the electron temperature T e tends to equilibrate 
rapidly with the vibrational temperature T v because of efficient vibrational tran- 
sitions by electron collisions. When T v becomes nearly equal to T e , the electron 
collisions cause no further change in vibrational populations, even though these 
temperatures may be far out of equilibrium with the translational temperature 
of the heavy particle T. After T e becomes equal to T,, heavy particle collisions 
promote T e and T v to equilibrate with the heavy particle translational temper- 
ature T (Translational- Vibrational energy transfer) 9 . Consequently, the above 
assumptions are validated at the initial stage of vibrational relaxation, where 
major transitions are confined to the relatively lower vibrational energy states, 
in gas mixtures by electron bombardment. 

Relaxation from Ground State (Case A) 


This is the case the flow immediately behind the shock wave where most 
molecules are in the ground state at t = 0. 

Lee 3,5 obtained the analytical solution for Eq. (6) assuming that the spatial 
variation of the second moment M is small. Then, the normalized number 
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density p v can be dictated as 

p v {v,t) = po[l - erf 


2 (tf 0 *) l/3 


( 10 ) 


where po = (n,/n» B ) v=0 is the normalized number density at v = 0 and K 0 = 
n t (M) v= o is the equivalent diffusivity at v = 0, and erf(z) is the error function 
defined by 


r /( z ) 


= 7 = fe-?dl 
V* Jo 


( 11 ) 


The value of po can be approximately assumed to be 1. 

The total vibrational energy E v per unit volume for harmonic oscillators is 
given by 

TO ^TJl 

E v = Y vhi/n v bPv ~hv vn v ePvdv (12) 

v=Q 

where h is Planck’s constant and v is the frequency of oscillators. Thus, 


E v (t) = - [1 - ~ T)> 

where A = hi/n\f(l — e~ a ), a = B v /T t and 6 = 2{Kot) 1 ^ . 

The time derivative of E v (t) can be dictated by 

dJ sr = AKMl - 


( 13 ) 


For small f, 


£?.(*) = A(iTo< H — ) 


( 14 ) 

( 15 ) 


For large t, 


dEJt) 

dt 


= AK 0 


E,(t) = A±( 1- 


a 2 


dE v {t) _ 2,4ir 0 
dt afcv/ir 


( 16 ) 

( 17 ) 

( 18 ) 


Relaxation from Excited State (Case B) 

In the flow region away from the shcok wave, the gas has already excited with 
the vibrational temperature T v . Therefore, to describe the relaxation process 
from this excited state, the excitation and de-excitation from/to all levels need 
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to be considered. Since the second moment of vibrational energy transfer M 
can not be considered constant in this case, the diffusion equation (6) must 
be solved numerically. Once the number density of vth excited state n„(t) is 
determined numerically, then, E v (t) and dE v (t)/dt are calculated. 

Since the diffusion equation (6) is the one-dimensional parabolic partial dif- 
ferential equation, the Crank-Nicolson scheme is used to solve the equation 
numerically. 


Rate Equation and Relaxation Time 


In the actual numerical simulation, it is quite convenient to use an appro- 
priate rate equation, which is compatible with the system of flow conservation 
equation, instead of using E v (t) or dE v (t)/dt obtained from the diffusion equa- 
tion. This rate equation corresponds to the source term for vibrational energy 
equation for a flowing gas. The resulting vibrational energy equation is used 
in a computaion of nonequilbrium multitemperature flow, along with the con- 
servation equations of mass, momentum, and energy. Therefore, a modified 
Landau-Teller type of rate equation, which approximate the rigorous solution of 
Eq. (6), is developed. In order to to derive the rate equation and the relaxation 
time, an assumption of Boltzmann distribution is made. That is, the vibra- 
tional energy levels of molecules are assumed to be populated by a Boltzmann 
distribution with an unique vibrational temperature T v . 

The rate equation for vibrational relaxation at any given instant may be 
defined in general as 

dE v (t) _ E v — E v , 
dt r« 

where f is a correction factor and is given by 


/ = ( 


Ey ~~ E v j 
e;-e v% o } 


( 20 ) 


where superscript and subscript 0 refer to the equilibrium value and the time 
i — 0. Arbitary parameter s in Eq. (20) has a positive value between 1 (Landau- 
Teller model) and 3.5 (diffusion model). 

Neglecting the correction factor in Eq. (19), which is only slightly different 
from unity, the vibrational relaxation time , r ei is given by 


e; - e v 

T ‘ dE v (t)/dt 


( 21 ) 


Relaxation from Ground State (Case A) 
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Since the molecules are all in the ground state at t — 0 , E v% q in Eq. (20) is 
zero and therefore, the solution of Eq. (19) is given by 


t 


In this case, r e can be approximately given as 

1 


(22) 


1 


jg 2 

T< — ^ dE 9 (t)/di) t ~° a 2 Ko (Q v /T e ) 2 Ko {0 v /T e ) 2 (l/2)n t Yl v »-o^o t v , {y ! ) 2 

(23) 

and the vibrational relaxation time r e is independent of the vibrational temper- 
ature T v . 

The vibrational relaxation time r e in Eq. (23) may be rewritten in the form 
p e r € (atm-sec) as 

_ ( 04 ) 

P,T ' ' (fc)’(i) 


Relaxation from Excited State (Case B) 


Based on the assumption that the vibrational levels are occupied by a Boltz- 
mann distribution with the vibrational temperature T v , the vibrational energy 
i?„ t o at t = 0 is given by 


E Vf0 — E v 


|t-o= E 0 — 


nhfhv 

e«*/r. - 1 


!?«(«sO) 


(25) 


The solution of Eq. (19) is 

E, = (E; - E 0 ){1 - [1 + (s - l)-] T * r } + E 0 ( 26 ) 

. The numerator of Eq. (21) is defined by 

e; -e v = E' v (T e ) - e 9 (t v ) = (27) 

The rate of change of vibrational energy d£ v (t)/dt is established by the excita- 
tion and de-excitation rates from/to all levels. For harmonic oscillators 


dE v (t) ' _ n nM hv | v — v | {( — —)&*', t> — ( ( 28 ) 

dt " ^ km tim 


v = l v 1 —0 


Thus, 


. _ * % fTl TTl 

*M1 = n e n M hv{\ - «*./*.)£ E I * - *' | - e~^ T ’k v<v ,} 

(29) 


v = 1 t/'=0 
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In compression flow (heating) case, the number density of i/th energy state 
is small compared with the one of rth state when v 1 > v, and the excita- 
tion rate coefficient k VfV * is much smaller than the de-excitation rate coeffi- 
cient k v i }V whenv > 1. Thus, the sum over v l in Eq. (29) may be chosen as 
v* = 0, l’ 2, • ■ • , v - 1. That is, the energy gains for all transitions to the state 
v comes from below, and the energy losses for all transitions from state v goes 
downward. Then, 


m v —l 


1M1 = n «n E( v - v ')( e ' 

di »=i «'=o , % 

(30) 

The vibrational relaxation time r t is dictated from Eq. (21) 


= — D 


where 


E: =1 E:r = l o(^ - v')(e-^/^k v> , v - 


( e 9./T. _ !)-i _ ( e «./T. _ !)-i 

D ~ 1 _ e -0./T. 


(31) 


(32) 


It must be noted that the relaxation time r« in this case is a function of both 
the vibrational temperature T v and the electron temperature T e . 

And 


p,r« = kT c D 


At T,,T„ >0,, 


1 


Er=iE:r = o(* - «')(c-’' tf ’ /T *^> - 


D ~ 


T„(T C -T.) 
61 


Thus, p.r* can be approximated by 
> iT e 


P'T, - 


T.(T.-T.) 


ET-iEritt*-'.. - + (0v/t v )(vK, v , - v>k v . tV )} 


(33) 


(34) 


(35) 


Results and Discussion 

The electronic ground state of the nitrogen molecule N 2 (X 1 E^ _ ) is consid- 
ered in the present study, which is the most predominant vibrationaily excited 
molecular species expected in the Earth-return APEV environment. Twelve 
initial vibrational states (v = 0 — 12) are included with 10 quanta change 
(Au = v r — v = ±10). Vibrational excitation and de-excitation rate coeffi- 
cients of N 2 , based on computational quantum chemistry 12 , are used in the 
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present calculation. These coefficients obtained by 5uo et. al. 13 provide only 
5 quanta change (At; = ±5) because of presumably precipitous dropping of the 
value for changes of 6 or more vibrational quanta. Hence, these rate coefficients 
are extrapolated to up to 10 quata change by using the exponential curve lit 
formula y = a • exp(-bx) in order to give sufficient convergence of the present 
calculation. The rotaional quantum number / = 50 is considered here as only 
the rate coefficients with this value are available for the present. 

Fig. 1 represents the e-V relaxation time in the form p t r t (atm-sec) for 
the case A (excitation only from the ground state) given by Eq. (24). As 
may be seen in Fig. 1, there is virtually no differnce between the result based 
on Lee’s semi-empirical rate coefficients 5 (dotted curve) and the one based on 
more rigorous coefficients of Huo et, al. 13 (solid curve). These theoretical 
rel ax ation times agree substantially with the result based on the experimental 
values of the cross sections of Schulz 7 , except at very low electron temperature 
(T f ~ 1000 K). It must be emphasized that Schulz’s experiment was conducted 
at room temperature, leading to the corresponding rotaional quantum number 
/ ~ 5. At low electron temperature, because of the small number density of the 
electrons with sufficient energy to cause resonant excitation, the contribution of 
nonresonant excitation (direct excitation) is dominant 5 . In both Lee 5 and Huo 
et. al. 13 calculations, this nonresonant direct excitation is neglected because of 
small effects on the excitation mechanism at moderate electron temperatures. 
Therefore, the experimental result is much more realistic at very low electron 
temperature. 

The e-V relaxation time for case B (excitation/de-excitation from the excited 
states) given by Eq. (33) is shown in Fig. 2. The relaxation time in this case is a 
function of both the electron temperature T e and the vibrational temperature T v 
[see Eq. (33)]. The vibrational temperature T v dictates the initial excited state 
with a Bolzmann distribution, and the molecules commence to equilibrate from 
T, to T«. As can be seen in Fig. 2, the relaxation times for the case B depend 
only weakly on the vibrational temperature T v . The relaxation times become 
slightly longer with the increse of the vibrational temperature T*, but the general 
tendency is exactly the same. The effects of multiple excitation/de-excitation 
from the excited states contribute to shortening the relaxation times noticeably 
in the entire temperature range as seen in Fig. 2. As can be deduced from 
the above-mentioned results, it may be concluded that the e-V relaxation times 
can be represented by an empirical formula which fits the Schulz’s experimental 
data at very low temperature and the equation (33) with a specific value of T v 
(say T v = 0.1 eV) at moderate and high temperatures. This formula will be a 
function of the electron temperature T e only, and applicable to the temperature 
range of 1000 K < T e < 50000 K. 

The second moment of vibrational transition rate coefficients M(v) is shown 
in Fig. 3 for the case where T e — 4000 K. The M(v) values are good indicator 
of the overall rate of transitions. M(v) increases rapidly from the ground state 


10 


value 3f(0) and remains nearly constant over the range of v = 1 to 12. This 
trend shows the strong effects of multiple excitation/de-excitation processes in 
the e-V vibrational transitions. This is also in contrast to the Landau-Teller 
model in which M(v) increases monotonically with the vibrational level v due 
to the fact that is proportional to v + 1 in the Landau-Teller model 11 . 

Figure 4 shows the normalized number density p v at selected times. This is a 
case wherein the molecules, initially in equilibrium at OAeV with a Boltzmann 
distribution, is heated suddenly by an electron heat bath having a constant 
temperature T e = 1.0 eV. As can be seen in Fig. 4, the population distribu- 
tion is quite different from a pure Boltzmann distribution in the early stage 
of relaxation process. Although the populations at v = 0 — 3 levels always 
maintain a Boltzmann distribution (a straight line in the v — p v plot), the pop- 
ulations at v > 4 considerably deviate from a Boltzmann distribution until at 
t = 5.48 x 10“ 11 sec. This underpopulation leads to suppress the equilibration 
process. Hence, the Boltzmann distribution assumption with unique T v may be 
questionable in the relaxation process from the excited states. 

The time variation of vibrational energy E v (t) given by the equation (12) 
is compared with the solution of rate equation (26) in Fig. 5 for the case of 
T,( f=0 ) = 0.4 eV ) T, = 1.0 eV and for n e = 10 l9 l/cm 3 . The abscissa in the fig- 
ure is the time normalized by the relaxation time r e [Eq. (31)]. The vibrational 
energy E v (t/r e ) is normalized by the equilibrium E*. As expected, the solid 
curve [Eq. (12)] reveals a different feature from the dashed curve [Eq. (26) with 
j = 3.5], due to the suppression phenomenon caused by a non-Boltzmann pop- 
ulation distribution. As evident in Fig. 4, the number densities of vibrational 
levels are underpopulated until t = 5.48 x 1G“ 11 sec, which is corresponding to 
t/r 9 = 5 in Fig. 5. Because of this under population, the vibrational energy is 
always smaller than the one predicted by a Boltzmann distribution during the 
early stage of relaxation. In the case A (relaxation from the ground state), the 
solution of rate equation (22) with s = 3.5 shows the diffusion characterise, 
and gives a good agreement with the solution of diffusion equation (13) in the 
temperature range of 0.1 eV < T e < 2.0 eV [see Ref. (5)]. However, in the case 
of B (relaxation from the excited states), the agreement is poor because of the 
profound effectes of underpopulated non-Boltzmann distribution in the early 
stage of relaxation. The development of more realistic rate equation is required. 


At present, effort is being made to obtain an empirical formula for the e-V 
relaxation time applicable to the temperature range of 1000 K < T t < 50000 K . 
The development of the rate equation, which represents the non-Boltzmann na- 
ture in the process of excitation/de-excitation from the excited states, is also 
underway. These results will be shown in the final paper along with the de- 
tailed discussion of the effects of non-Boltzmann distribution on the vibrational 
equilibration process. 
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